The torque transfer coefficient in DNA under torsional stress 
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In recent years, significant progress in understanding the properties of supercoiled DNA has 
been obtained due to nanotechniques that made stretching and twisting of single molecules possi- 
ble. Quantitative interpretation of such experiments requires accurate knowledge of torques inside 
manipulated DNA. This paper argues that it is not possible to transfer the entire magnitudes of 
external torques to the twisting stress of the double helix, and that a reducing torque transfer coeffi- 
cient (TTC<1) should always be assumed. This assertion agrees with simple physical intuition and 
is supported by the results of all-atom molecular dynamics (MD) simulations. According to MD, 
the TTCs around 0.8 are observed in nearly optimal conditions. Reaching higher values requires 
special efforts and it should be difficult in practice. The TTC can be partially responsible for the 
persistent discrepancies between the twisting rigidity of DNA measured by different methods. 

PACS numbers: 87.14.gk 87.15.H- 87.15. ap 87.15. ak 

The twisting stiffness is a unique property of DNA 
that makes possible its supercoiling, which is essential 
for genome compaction and regulation [![. At present, 
the mechanics of supercoiled DNA is frequently studied 
by using nanotechniques that offer a means to stretch 
and twist single molecules 0-0] • Due to its remark- 
able conceptual simplicity this method has gained broad 
recognition. The effects of external twisting depend upon 
the torsional rigidity of DNA and the values of applied 
torques. Accurate knowledge of these quantities is neces- 
sary for interpretation of experimental data. The recent 
methodological advances made twisting of single DNA 
with external torques of known magnitude [H, la] possi- 
ble. However, one question has never been asked: Is it 
evident that the twisting stress created inside DNA al- 
ways corresponds to the external torque measured from 
outside? 

A simple example shown in Fig. Q] demonstrates that 
the answer to the above question can be negative. Con- 
sider the left panel. An elastic rod is fixed on a solid base, 
with a constant torque r applied to its top, which causes 
a deviation of the equilibrium twist angle A$ r = $ r — $o- 
For simplicity, assume that the bending rigidity is very 
high and bending deformations can be neglected. If the 
torque axis is deviated by angle 9 only the vertical torque 
projection works and the twisting stress inside the cylin- 
der as well as A$ r decrease as cos#. The above effect 
persists with the addition of other degrees of freedom. 
Fig. Q] illustrates this for Brownian dynamics (BD) of 
a discrete wormlike rod (WLR) model where bending, 
twisting and translational motions of DNA are repro- 
duced in agreement with experiments Q. In the har- 
monic approximation the torsional part of the free energy 
is 
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FIG. 1: (Color online) (Left panel) An elastic rod is fixed on 
a solid base. A twisting torque is applied at the top with 
the torque axis inclined by angle 9 form the vertical axis. 
(Right panel) A discrete WLR model of a DNA fragment of 
three base pairs. The external torques are applied to terminal 
bases under angle 9 with respect to the chain axis. (Middle 
panel) Twisting caused by external torques in BD of a 14- 
mer DNA fragment with elasticity corresponding to MD data 
0. The torques were applied with 9 = 0° (black dots), 30° 
(red squares), and 60° (blue circles). The black dashed line 
displays the theoretical dependence for the input value of It- 



and T is the temperature. The twisting rigidity of DNA 
is characterized by parameter l t called the torsional per- 
sistence length. In this case 



A<i% 



tL 
kfft 



(2) 



and the probability distribution of torsional fluctuations 
is Gaussian 
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where L is the DNA length, k is Boltzmann's constant, 



The right panel of Fig. [T] shows a coarse-grained WLR 
representation of a DNA trimer. Each base pair is mod- 
eled by a rigid composite bead of four particles O, X, 
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Y, and Z that define a local Cartesian frame. The OZ 
axis is constrained to pass through the center of the 
following bead in the chain. Other degrees of freedom 
are restrained by harmonic nearest-neighbor potentials 
adjusted to provide the desired macroscopic rigidity for 
bending and twisting The twisting torques are ap- 
plied to terminal composite beads by using the earlier 
algorithm with a small modification to deviate the 
torque vector by a fixed angle 8 with respect to the local 
OZ axes. The external torque at one end is compensated 
by reactions at the opposite end so that the integral ex- 
ternal force and torque on the molecule are always zero. 
This can be viewed as if a demon in the moving global 
frame twists one end of the molecule and holds the oppo- 
site end so that the overall translation and rotation are 
not perturbed. The results of BD simulations of a 14- 
mer DNA fragment are shown in the middle panel. As 
expected, the effective torque decreases approximately as 
cos 9. The analytical linear dependence corresponding to 
the true It value is reproduced in BD only with 0=0. 
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FIG. 2: (Color online) The normalized probability distribu- 
tions of twisting fluctuations in DNA obtained by all-atom 
MD simulations 0. The $ angles are measured directly be- 
tween the terminal base pair frames; therefore, $ = for one 
helical turn. The CG dodecamer was considered in relaxed 
state (r = 0) and under steady torsional stress induced by 
external torques of r = -80 and +60 pN nm, respectively. 
The dashed lines show the corresponding Gaussian distribu- 
tions [Eq. (|3|] for the computed values of $ T and l t . The 
solid lines display the analytical distributions for the best-fit 
anharmonic torsional potential 

It might seem that the above example is irrelevant be- 
cause in experiments the external torques are always ap- 
plied along the stretched DNA double helix. The prob- 
lem is that local twisting axes in real DNA may be in- 
clined with respect to the helix direction. This issue sur- 
faced unexpectedly in the recent study of elastic prop- 
erties of strained DNA by all-atom molecular dynam- 
ics (MD) 0, E3]. In these computations steady twist- 
ing torques were applied to two short DNA fragments 
that I denote here as TA and CG (because of the AT- 
and GC-alternating sequences, respectively). The clas- 
sical all-atom MD produce the probability distribution 
of twisting P$ . The three distributions shown in Fig. m 
were obtained for dodecamer CG (one helical turn) with 
contrasting values of the external torque. From these 
data, the It value can be extracted by using either Eq. 



(0) or and the resulting two values should be iden- 
tical. However, according to Eq. ([3]) the variance of 
the distribution should be constant whereas in Fig. [2] it 
evidently changes with r; that is, the harmonic approx- 
imation breaks down. Therefore, a broader issue should 
be considered, namely, whether or not these MD results 
agree with statistical mechanics of an anharmonic elastic 
rod. To answer this question we look for an anharmonic 
potential that would consistently describe the MD data. 
Assuming that <5> T = for zero load the free energy is 

f/ T ($) = -t$ + c7($), 

where J7($) is the torsional potential of bare DNA, there- 
fore, 
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According to Fig. [2] all the probability distributions are 
nearly Gaussian; that is, in the vicinity of $ T the har- 
monic approximation is valid and one can use the Taylor 
series to obtain 
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Let us look for a polynomial form of £/($). We choose 
£/(0)=0; in this case the first term is quadratic. For i t (r) 
not to be constant the polynomial should be of the third 
degree or higher, for instance, 
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In this case Eq. ((4j and dU give 

1 o 1 , 
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MD provides $ r and It measured for several values of r. 
Coefficients qi are found by the least squares minimiza- 
tion of the discrepancy between the left- and right-hand 
parts of Eq. ([7]) and (|SJ). The results of this fitting are 
shown in Fig. |3]by dashed red lines. It is seen that they 
significantly deviate from the MD points. The deviations 
are clearly systematic and they go far beyond the error 
ranges. This situation does not change with increased 
degree of the polynomial. Note that in all plots the de- 
viations grow with |t|. It turned out that the quality of 
fitting is radically improved if a torque transfer coefficient 
(TTC) q is added to Eq. Q as follows: 



1 9 1 1 

q T = qi <P T + q 2 -^ T + q 3 -K, 



(9) 



with qo optimized together with other qi. In this case, the 
agreement within the range of statistical errors is easily 
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reached for polynomials of the fourth degree and higher. 
Interestingly, the optimized TTC values are similar for 
CG and TA (0.81 and 0.79, respectively) and they vary 
by only a few percent when the degree of the polynomial 
is increased up to six. The quality of fitting is illustrated 
by the solid traces in Fig. [5]and[3] The agreement is quite 
good, notably, small deviations of P$ from Gaussians in 
Fig. [2]are well reproduced. These subtleties have low sta- 
tistical weight and actually do not affect the values of $ r 
and It involved in the fitting. One can conclude that the 
MD results qualitatively and quantitatively agree with 
statistical mechanics, but only if the occurrence of the 
TTC qo<l in Eq. (|9]) can be reasonably explained. 
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FIG. 3: (Color online) Fitting of analytical models to the 
results of all-atom MD simulations of two DNA fragments 
(TA and CG) [7|. Error bars show statistical errors. In the 
upper two panels the symbol size corresponds to maximal 
errors. The dashed red lines show the results of fitting for a 
polynomial torsional potential without a TTC applied. The 
solid black lines show best fit dependences for a fourth order 
polynomial with optimized TTC. 

The difference between Eq. (J7J, and © indicates that 
good fitting requires that the values in Eq. (J7]) and flSJ) 
differ by a constant factor. Therefore, the scaling might 
be transferred from Eq. ^ to Eq. (JSJ , which would give 
equally good fitting. At first sight, this seems reasonable 
because the scaling can be interpreted as a correction 
of the uncertainty in the L/T value. Actually, however, 
this factor is used in computing l t from MD data and 
a correcting scaling should be applied simultaneously to 
both sides of Eq. ©. Therefore, the occurrence of (ft can 
be reasonably explained only in Eq. © by assuming that 
a fraction of external torque applied to DNA is somehow 
lost. 

The possible origin of incomplete torque transfer was 
shown in Fig. [TJ above. The DNA can respond to exter- 



nal torques by either twisting or bending or both. The 
torque fraction that causes twisting varies with angle 9 
whereas the remaining part is largely canceled by auxil- 
iary reaction forces. Such forces are always necessary. In 
the left panel of Fig. [TJand in single molecule experiments 
the reaction forces maintain DNA straight and fixed in 
place. In computer simulations, they are used to zero the 
integral external force and torque that otherwise would 
accelerate the overall translation and rotation Q. All 
these reactions do not affect torsional fluctuations, but 
effectively damp all other deformations, which is where a 
part of the external torque goes. 

The orientation of the twisting axis in real DNA is 
not known. The external torques in MD are applied to 
terminal bases, with the torque axes orthogonal to the 
base plane Q. On closer inspection, however, one notices 
that this choice may not be optimal. There are other 
distinguished directions, for instance, the helical axis or 
the perpendicular to the base-pair plane. However, these 
directions are well defined and coincident only in the ideal 
B-DNA model. In reality, bases are not perpendicular to 
the helix direction and paired bases are not co-planar, 
therefore, the definition of the base pair plane as well as 
the helical axis is a matter of convention. Moreover, local 
twisting occurs due to many rotatable covalent bonds and 
it is possible that the transfer of the torsional stress to 
DNA depends upon specific orientations of these bonds 
near the point where the torque is applied. 

To get an idea of the situation in real DNA, I tried to 
vary the 9 angle in MD similarly to the BD simulations 
presented in Fig. [TJ The protocols of MD simulations 
were the same as before 0, [l(| ■ The external torque is 
applied by using a Cartesian frame rigidly attached to 
the base [9JJ. The axes of this frame approximately cor- 
respond to the standard convention (see Fig. [3}, an d by 
default the torque is directed along the Z axis. Unlike the 
WLR model, the all-atom DNA does not have cylindri- 
cal symmetry, therefore with 9 ^ 0, we should also check 
different torque azimuths. The natural qualitatively dis- 
tinct directions of deviation correspond to XZ and Y Z 
planes. With 6 ^ in the XZ plane, the torque vector 
is inclined toward the major {9 < 0) or minor (9 > 0) 
groove. In the YZ plane, the torque can be inclined to- 
ward the backbone of the same (9 < 0) or the opposite 
(9 > 0) base. The magnitude of the external torque was 
fixed at ±40 pN nm. A relatively large value is required 
to obtain variations beyond the range of statistical er- 
rors. Earlier studies showed that, in MD, the properties 
of short DNA fragments change smoothly in this range 
of torques 

The results shown in Fig. [4] confirm that the effect 
is also significant in the all-atom MD. Computed points 
are compared with theoretical dependences of the form 
Scos(9 — 9 S ), with 5 and 9 S adjusted manually. Angle 
9 S specifies the softest local twisting axis. Parameter 8 
affects simultaneously the vertical displacement and the 
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FIG. 4: (Color online) The upper panel displays the orien- 
tation of t he p rincipal local axes according to the standard 
convention The Z axis is directed along the double he- 
lix. The X axis looks towards the major groove. The Y axis 
is directed along the base pair so that the three axes form a 
right-hand triple. The torque vector is inclined by a variable 
angle 9 with respect to the Z axis. The lower panel displays 
the results of all-atom MD simulations of DNA subjected to 
external torques r = ±40 pN nm applied to terminal bases. 
Red open circles show the results for the torque vector inclined 
in the XZ plane (toward the grooves), with 9 > correspond- 
ing to deviations toward the minor groove. Similar data for 
the YZ plane (deviations toward the backbones) are shown 
by blue open squares. In this case, the positive 9 values corre- 
spond to deviations toward the opposite base. The solid lines 
display the plots of Scos(9 — 9 a ), with 8 and 9 S fit manually 
to match the computed points. 



inclination of the curves, i.e., the fitting is constrained. 
The agreement is relatively good for twisting and less so 
for untwisting. In both planes the softest axes are found 
close to, but not exactly perpendicular to the base. For 
deviations in the Y Z plane the optimal 9 S value is about 
-5° for both signs of the torque. In the XZ plane the 
softest twisting direction is observed at 9 S = +10° and 
+20° for twisting and untwisting, respectively. There- 
fore, the all-azimuth softest axis may be inclined by 20- 
25° toward the minor groove, which would correspond to 
a TTC value around 0.9. Taking into account the de- 
tails of the original algorithm the TTC can be reduced 
to 0.87 12|. The remaining discrepancy with the best 
fit value qo « 0.8 is not negligible. Moreover, when the 



twisting is measured by summing the local twist angle 
along the chain the optimized TTC values decrease 
to 0.72 and 0.76 for TA and CG fragments, respectively. 
This residual discrepancy can be attributed to confor- 
mational fluctuations. In the ideal WLR model and BD 
simulations shown in Fig. [TJ the fixed angle 9 is mea- 
sured between the torque and the axis of local twisting. 



In contrast, in Fig. @] and MD simulations, the exter- 
nal torque is applied with 9 fixed with respect to the 
base plane whereas the local twisting axis may fluctuate 
with the DNA structure. As a result, even in an optimal 
orientation, the torque is not strictly collinear with the 
twisting axis, and this should further reduce the TTC. 

The foregoing results prove that it is difficult, if ever 
possible, to transfer the desired value of external torque 
to real double helical DNA. This conclusion should 
also be valid for nanomanipulations with isolated DNA 
molecules. The samples used in experiments commonly 
involve a few kilobase pairs of nearly random DNA. Ear- 
lier studies indicate that the WLR model is valid for any 
integral number of helical turns starting from one [illl3j: 
therefore, one helical turn used in MD already involves all 
essential properties and the conclusions obtained should 
be invariant to length scaling. In experiments the double 
helix is attached to macroscopic beads and extended by a 
stretching force. The torques are applied to the beads in 
different ways and transferred to DNA through a num- 
ber of point contacts. However, any multi-point interac- 
tion can be decomposed to individual components, each 
characterized by a specific TTC. The quasi-randomness 
of natural DNA sequences and the possible sequence de- 
pendence cannot smooth out the effect because the TTC 
never exceeds 1. Therefore, the induced twisting is al- 
ways downgraded. 

The TTC can explain some discrepancies in the re- 
cent experimental measurements of the torsional rigidity 
of DNA. The value of It (or C ~ kTl t ) has been mea- 
sured since the 1970s [3] because it plays a fundamental 
role in all manifestations of DNA supercoiling. Before the 
recent single-molecule studies, C ss 310 pN nm 2 was con- 
sidered as a consensus value for random sequence DNA. 
This value reasonably agrees with some indirect estimates 



that used single-molecule data [15|, [16| , but it is at least 
20-30% lower than in direct measurements using forced 
DNA twisting by calibrated external torques 0, @| . Sev- 
eral pitfalls possibly responsible for such a difference have 
been discussed in the literature 0, [l?}- Assuming the 
TTC values around 0.8, which may slightly vary between 
different experimental installations, all these data can be 
accounted for without additional assumptions. Only one 
observation cannot be explained, namely, the very high 
apparent torsional rigidity (C = 440 pN nm 2 ) reported 
for equilibrium thermal torsional fluctuations of single 
stretched DNA [Sj]. I believe that some other factors 
could affect this result and it requires further considera- 
tion. 

In summary, the present paper shows that because of 
the structural details of the double helix, it is not possi- 
ble to transfer the entire magnitudes of external torques 
to the internal twisting stress of DNA, and that a reduc- 
ing torque transfer coefficient should always be assumed. 
This assertion agrees with simple physical intuition and 
it is supported by the results of all-atom MD simulations. 
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Theoretical estimates indicate that the TTCs around 0.8 
should be observed for torques nearly parallel to the dou- 
ble helix. Reaching higher values requires special efforts 
that should take into account the details of the DNA 
structure and it should be very difficult for experiments 
with long molecules. The TTC explains some controver- 
sies in the experimental values of DNA twisting rigidity 
measured by different methods. 

The author is grateful to Alex Vologodskii for useful 
discussions and comments to the paper. I also thank 
Zev Bryant for correspondence and valuable clarifications 
concerning the rotor bead tracking method. 
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